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D. H. Mcintosh 2 , H.-W. Rix 3 , M. J. Rieke 2 , C. B. Foltz 4 

ABSTRACT 

We have observed a sample of 22 luminous quasars, in the range 2.0 <^ z <^ 2.5, 
at 1.6/im with the near-infrared (NIR) spectrograph FSPEC on the Multiple Mirror 
Telescope. Our sample contains 13 radio-loud and 9 radio-quiet objects. We have 
measured the systemic redshifts z sys directly from the strong [O III]A5007 line 
emitted from the narrow-line-region. From the same spectra, we have found that the 
non-resonance broad H/3 lines have a systematic mean redward shift of 520 ± 80 km s _1 
with respect to systemic. Such a shift was not found in our identical analysis of the 
low-redshift sample of Boroson & Green. The amplitude of this redshift is comparable 
to half the expected gravitational redshift and transverse Doppler effects, and is 
consistent with a correlation between redshift differences and quasar luminosity. From 
data in the literature, we confirm that the high-ionization rest-frame ultraviolet broad 
lines are blueshifted ~ 550-1050 km s^ 1 from systemic, and that these velocity shifts 
systematically increase with ionization potential. Our results allow us to quantify the 
known bias in estimating the ionizing flux from the inter-galactic-medium J^ GM via 
the Proximity Effect. Using redshift measurements commonly determined from strong 
broad line species, like Lya or C IVA1549, results in an over-estimation of J^ GM by 
factors of ~ 1.9 — 2.3. Similarly, corresponding lower limits on the density of baryons 
rij, will be over-estimated by factors of ~ 1.4 — 1.5. However, the low-ionization 
Mg IIA2798 broad line is within ~ 50 km s _1 of systemic, and thus would be the 
line of choice for determining the true redshift of 1.0 < z < 2.2 quasars without NIR 
spectroscopy, and z > 3.1 objects using NIR spectroscopy. 



Subject headings: cosmology: miscellaneous — infrared: general — quasars: emission 
lines — quasars: general — relativity 
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1. Introduction 



Defining the center-of-mass (COM = systemic) rest-frame, thus z sys , of quasars (QSOs) and 
of their host galaxies is important for a number of reasons. First, it is crucial for estimating 
the ionizing flux from the inter-galactic-medium J* GM via the Proximity Effect flBajtlik, Duncan 
fc Ostriker 1988 ). To estimate the redshift range (z S ys — -ZLyo) over which Lya absorbers are 



underpopulated, we need to know z sys accurately. An error of 1000 km s 1 in z sys leads to a 



factor of ~ 2.6 difference in the inferred J^ GM ( Espey 1993 ). In turn, precise estimates of J* M 
are needed to understand the ionization balance of the IGM and the structure of the Lya Forest. 
Furthermore, accurate knowledge of the systemic frame is important to understanding the 'central 
engine' and the broad-line-region (BLR) physics. Inflow, outflow and radiative transfer effects may 
lead to offsets between the emission line centroids and z sys . The value of z sys also plays a critical 



role in understanding the kinematics of 'associated' (-z a bs ~ z cm) absorption line systems QFoltz 
et al. 1986 ) and the nature of the 'z a b s > -Zem' metal line systems ( paskell 1982 ). In both cases, 



determination of z sys is necessary to understand whether the absorbing gas is in- or out-flowing, 
or at rest with respect to the QSO. 

As Tytler & Fan (1992) first showed, the QSO redshift as measured by different broad emission 
lines systematically depends on ionization potential. In general, it has been well documented 
(Gaskcll 1982; Espey et al. 1989| ; Brotherton et al. 1994; and references therein) that broad 
emission lines from highly ionized species (HILs; e.g. C IVA1549, C III]A1909, N VA1240), as well 
as Lya, are systematically blueshifted 500-1500 km s _1 with respect to broad emission lines from 
lower ionization species (LILs; e.g. Mg IIA2798, O IA1305), and the permitted H I Balmer series. 
In contrast, redshifts from narrow emission lines (e.g. [O III]A5007) in low redshift (z < 1) QSOs 
and Seyfert galaxies are observed to be within 100 km s _1 of the broad Mg II and Balmer lines 
QGaskell 1982| ; |Tytler fc Fan 1992| ; Poroson fc Green 1992| ; |Laor et al. 1995| ; and |Corbin fe Boroson 
1996 ). Furthermore, in local AGN, redshifts from narrow forbidden lines have shown agreement 



to <^ 100 km s _1 of the accepted systemic frame determined by stellar absorption features and H I 
21 cm emission in the host galaxies ( Gaskcll 1982j ; Vrtilek fc Carleton 1985 , Hutchings, Gower & 
Price 1987 ), and are within ±200 km s" 1 of the CO emission (|Alloin et al. 1992Q . Since the narrow 
emission originates from a large volume of gas centered on the QSO (i.e. the narrow-line-region 
with tnlr ~ 1 kpc), Penston (1977) first proposed that it is likely that the narrow forbidden lines 
give the COM redshift. 

At higher redshifts (z > 1.5) the strong narrow lines can only be detected in the NIR. Early 
NIR spectroscopy of a single bright z = 2.09 QSO showed Ha, H/3 and narrow [O III] emission 
had similar redshifts to the LILs, supporting the low redshift studies that claimed Balmer lines 
provided a reliable measure of z sys flCarswell et al. 1991 ). However, Nishihara et al. (1997) took 
NIR spectra of four z > 1.5 QSOs and found the broad Ha line to be redshifted on average by 
430 ± 170 km s^ 1 from the systemic [O III] narrow line. 



From our total sample of 32 bright (V < 18.0 mag), 2.0 £ z £ 2.5 QSOs with NIR if -band 
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spectra (see Mcintosh, Rieke, Rix, Foltz Sz Weymann 1999| ; hereafter M99), we selected the 
set of 22 objects with detected [0 III]A5007 emission ([O III] EW > 5A, except Q0842+345, 
Q1011+091, and Q1158-187 for which the low S/N prohibited reliable line identification). This 
nearly unprecedented detection of strong emission lines from the NLR of distant, high luminosity 
QSOs has provided the direct measurement of z sys via the rest-frame [O III]A5007 line center. 
Consequently, this unique data set allows direct comparisons between the COM frame and a 
variety of resonance and recombination broad line centers in this sample. In addition, each 
spectrum contains the broad H/3 line, thus, we are able to directly compare the NLR redshift with 
a redshift obtained from BLR emission in the same spectrum. 



2. Emission Line Redshift Measurements 

To measure z sys for each QSO, we fit roughly the upper half of the [O III]A5007 line with a 
Gaussian, except for three objects (Q1228+077, Q1331+170 and Q1416+091) for which the peak 
or centroid of the [O III] line was used due to the poor S/N of these spectra. This fitting procedure 
provided line center measurements with much smaller uncertainties compared with the broad H/3 
measurements (see below). We estimated the uncertainties to be a ~ ±50 km s _1 . The measured 
systemic redshifts are presented in column (2) of Table 1. 

To determine the QSO broad H/3 redshift z em , we shifted each spectrum to its systemic 
frame and fit it with a set of templates representing [O III], H/3, Fe II and continuum emission 
features using a \ 2 optimization routine (as described in M99). All components except H/3 were 
subtracted from the data spectrum resulting in an H/3 line spectrum. The central wavelength of 
the broad component was obtained by fitting a single Gaussian to all portions of the line profile 
with relative flux < | of the maximum (see M99). The broad H/3 line centers measured in the 
systemic frame, with their ±lcr uncertainties, are given in column (3) of Table 1. In addition, the 
emission redshifts derived from broad H/3 are listed in column (4). 
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Table 1. Observed Redshift Parameters 



QSO z sys Broad H/3 



Actr Zcm Au 









(A) 








(km 


s 




BAL0043+008 


2 


.146 


4866 


± 


2 


2 


149 


310 


± 


120 


Q0049+014 


2 


.307 


4862 


± 


4 


2 


308 


60 


± 


310 


Q0109+022 


2 


.351 


4868 


± 


6 


2 


356 


430 


± 


370 


Q0123+257 


2. 


370 


4868 


± 


3 


2 


.375 


430 


± 


190 


Q0153+744 


2 


341 


4871 


± 


2 


2 


.348 


620 


± 


190 


Q0226-038 


2 


.073 


4858 


± 


5 


2 


071 


-190 


± 


310 


Q0421+019 


2. 


.056 


4873 


± 


3 


2 


064 


740 


± 


190 


Q0424-131 


2. 


168 


4866 


± 


3 


2 


.171 


310 


± 


190 


Q0552+398 


2. 


.363 


4867 


± 


3 


2 


367 


370 


± 


190 


HE1104-181 


2. 


318 


4870 


± 


3 


2 


.324 


560 


± 


190 


Q1148-001 


1 


980 a 


















Q1222+228 


2 


.058 


4858 


± 


3 


2 


.056 


-190 


± 


190 


Q1228+077 


2 


.389 


4860 


± 


3 


2 


388 


-60 


± 


190 


PG1247+267 


2. 


042 


4873 


± 


1 


2 


.050 


740 


± 


60 


Q1331+170 


2. 


.097 


4861 


± 


4 


2 


097 





± 


250 


Q1416+091 


2. 


.017 


4851 


± 


8 


2 


Oil 


-620 


± 


490 


Q1435+638 


2. 


066 


4875 


± 


4 


2 


075 


860 


± 


250 


Q 1448-232 


2 


.220 


4876 


± 


2 


2 


.230 


930 


± 


120 


Q1704+710 


2 


.010 


4879 


± 


5 


2 


.021 


1110 


± 


310 


BAL2212-179 


2. 


.228 


4864 


± 


5 


2 


.230 


190 


± 


310 


Q2251+244 


2. 


359 


4852 


± 


4 


2 


.353 


-560 


± 


250 


Q2310+385 


2 


.181 


4869 


± 


4 


2 


186 


490 


± 


250 



a Redshift of this QSO is such that H/3 line center is 
blucward of //-band window. 



- 5 - 



3. Velocity Shifts 



Comparing the z sys and the broad H/3 z em measurements for each QSO (see Table 1) reveals 
a velocity shift Av relative to the systemic frame. A Av > corresponds to a redward shift of H/3 
with respect to systemic, while a negative At; represents a blueshift. The velocity shifts for broad 
H/3 are tabulated in column (5) of Table 1, where the uncertainties are determined by the broad 
HP component. 

This is the first sizeable sample of high luminosity QSOs to permit a statistical comparison 
of the velocity shifts of various UV rest-frame broad line redshifts relative to systemic. Therefore, 
we compiled the published redshifts for each object in Table 2, and calculated At; for each line 
species: 



where zn ne is the published redshift. The distribution of velocity shifts for each UV line species is 
plotted separately in Figure 1. All blueshifts larger than 2250 km s _1 are due to a single radio-loud 
object (Q2251+244) with strong [O III] emission (see M99, Fig. 1). The z sys measurements are 
accurate to ~ 50 km s _1 on average, while the majority of the published redshifts are known to 
±0.001 which corresponds to roughly ±100 km s _1 . 




(1) 
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Table 2. Published Emission Line Redshifts 



QSO 


Lya 
A1216 


Refs. a 


nV 

A1240 


Refs. a Si 


IV/O IV] b 
A1400 


Refs. a 


cIV 

A 1549 


Refs. a 


cIII] 

A1909 


Refs. a 


Mg II 

A2798 


Rcfs. a 


BAL0043+008 














2 


141 


1,2 


2.143 


2,3 


2.143 


3 C 


Q0049+014 


2.31 


4 










2 


285 


5 


2.30 


6 d 






Q0109+022 


2.35 


4 










2 


348 


5 


2.34 


6 d 






Q0123+257 


2.360 


7,8 d 










2 


356 


7,8 d 


2.36 


6 d 






Q0153+744 


2.342 


9 


2.340 


9 


2.340 


9 


2 


339 


9 


2.339 


9 


2.349 


9 


Q0226-038 


2.069 


10-12 


2.058 


12 


2.057 


10,12,13 


2 


064 


10-14 


2.064 


14 


2.073 


14 


Q0421+019 


2.048 


12,15 










2 


051 


15 


2.053 


12,14 


2.056 


14 


Q0424-131 


2.160 


13,16 






2.168 


13 


2 


164 


13,14,16 


2.164 


14,16 


2.166 


14,16 


Q0552+398 


2.362 


4,17 
























HE1104-181 


2.315 


18 










2 


312 


18-20 










Q1148-001 


1.979 


11,12 






1.975 


10,11,12 


1 


979 


10-14,21 


1.979 


10,14,21 


1.978 


14,21 


Q1222+228 


2.046 


11,12,16 


2.049 


12 


2.046 


12,13 


2 


042 


11-13,16 


2.053 


16 


2.056 


16 


Q1228+077 


2.388 


4,22,23 










2 


389 


24 


2.392 


12 






PG1247+267 


2.040 


11-13 


2.042 


12 


2.040 


12,13 


2 


039 


11-13 


2.043 


25 d 






Q1331+170 


2.084 


11,12,16,26 


2.079 


12,26 


2.086 


12,13 


2 


079 


11-13,16,26 


2.086 


16,26 


2.095 


16,26 


Q1416+091 




















2.012 


14 


2.018 


14 


Q1435+638 










2.072 


13 


2 


066 


13,14 


2.062 


14,27 


2.061 


14 


Q1448-232 


2.216 


16,28 










2 


216 


16,28,29 


2.214 


16,29 


2.223 


16 


Q1704+710 




















2.011 


14 


2.019 


14 


BAL2212-179 


2.220 


30° 


2.205 


30° 


2.218 


30° 


2 


209 


30° 










Q2251+244 


2.330 


7,10,31 






2.328 


10,12 


2 


323 


7,10,12,31 


2.318 


12,31 






Q2310+385 














2 


175 


32 










a Multiplc references indicate that the rcdshift is an average. 
b A1400 blend of the Si IV doublet and the O IV multiplet. 
c Rcdshift extrapolated from an emission line profile comparison. 



Rcdshift measured directly from published spectrum. 
°Redshift measured directly from digital spectrum. 



References. — (1) Osmcr et al. 1994; (2) Turnshek et al. 1980; (3) Hartig & Baldwin 1986; (4) Wolfe et al. 1986; (5) Schneider et al. 1994; (6) Pci et al. 1991; 
(7) Schmidt 1975; (8) Barlow & Sargent 1997; (9) Lawrence et al. 1996; (10) Gaskell 1982; (11) Young et al. 1982; (12) Tytlcr & Fan 1992; (13) Sargent et al. 
1988; (14) Stcidel & Sargent 1991; (15) Schmidt 1977; (16) Espey et al. 1989; (17) Wills & Wills 1976; (18) Smctte et al. 1995; (19) Wisotzki et al. 1993; (20) 
Reimers et al. 1996; (21) Aldcroft et al. 1994; (22) Vaucher & Wcedman 1980; (23) Robertson & Shaver 1983; (24) Sramek & Weedman 1978; (25) Green et al. 
1980; (26) Carswell et al. 1991; (27) Laor et al. 1995; (28) Jian-sheng et al. 1984; (29) Ulrich 1989; (30) Morris et al. 1991; (31) Barthcl et al. 1990; (32) B. J. 
Wills 1997, private communication 
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4. Results and Discussion 

From the above analysis, we find: 



that the H/3 broad component emission has a weighted mean redshift of 520 ± 80 km s 1 
relative to z sys , with a sample variance of 360 km s _1 . 

that on average the Mg IIA2798 based redshift is within ~ 50 km s _1 of systemic, confirming 



that LILs should provide reliable z sys estimates ( Carswell et al. 1991 ). The rest wavelength 
and strength of Mg II emission make it the line of choice for determining z sys of distant 
(1.0 < z < 2.2) QSOs without NIR spectroscopy, or even more remote [z p> 3.1) objects 
using NIR spectroscopy. 

3. confirmation of the following effects: (i) The average blueshift (~ 600 - 1100 km s _1 ) of the 
HILs relative to the LIL Mg II. (ii) The correlation between average blueward velocity shift 
and ionization potential for each resonance broad line species. 

4. no correlation between QSO radio type (i.e. radio-loud versus radio-quiet) and any of the 
line species velocity shifts. 



4.1. Redshifted Balmer Lines 

The mean redshift of broad H/3 is similar to that found for Ha in a small sample (4) of 
luminous z ~ 2 QSOs ( [Nishihara et al. 1997 ). Yet at lower redshifts (z < 1), Balmer emission 



gives the same redshift to within 100 km s~ x of narrow forbidden emission. Given that z ~ 2 
QSOs are more luminous than their present-day counterparts, perhaps this observed trend of 
increased Balmer redshift with increased z SJS represents a luminosity dependence. In fact, we 
find a correlation between H/3 redshift and QSO rest-frame y-band luminosity (from M99), such 
that the ~ 100 times less luminous sample of Boroson & Green (1992) exhibits a mean broad H/3 
redshift of < 100 km s _1 (see Figure 2). We conclude that the use of Balmer emission lines to 
measure z sys may not be prudent for more distant, and thus more luminous QSOs. 



The magnitude of the redward Balmer shift may be partially due to relativistic effects flNetzer 



1977). Recent studies (e.g. Corbin fc Boroson 1996j ) have suggested a connection between the 



luminosity dependence of QSO broad line redward asymmetries and gravitational redshifts due 
to virialized motions near super-massive (~ 10 9 -10 10 Mq) black holes. Assuming BLR emission 
from a gravitationally bound (virialized), spherical ensemble of clouds at r ~ 0.1 pc allows the 
conversion from the observed mean broad line width to the mean line-of-sight (LOS) velocity and 
the true average cloud velocity: (f fwhm) = 2(wlos) = ^(^cioud)- We measured the weighted 
mean line width (dfwhm) = 8750 ± 570 km s _1 , with a sample variance of 2550 km s _1 , from the 
H/3 broad component fitting (M99). The gravitational redshifting of emission from clouds deep 
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within the potential well of the super-massive black hole is Augr — cl g ud — 190 ± 20 km s . In 
addition, emission from clouds moving perpendicular to our LOS experiences a transverse Doppler 
shift of Atjtds ~ ^ ~ ~ 60 ± 10 km s _1 . 

The combination of these two relativistic effects can account for roughly half of the 
520 ± 80 km s™ 1 Balmer redshift. Although intriguing, it is unclear why this emission would not 
also exhibit Doppler beaming effects (e.g. [Laor 1991] ) and why the HILs, which reverberation 
mapping studies have placed at the same distance, have blueshifts. 



4.2. Cosmological Implications of the z sys Measurements 

HILs, like Lya or C IV, are not directly measuring z sys . Specifically, for this sample we find 
Lya and C IV to be blueshifted on average 550 and 860 km s _1 , respectively, relative to systemic. 
Therefore, using redshifts acquired from these lines will result in the over-estimation of J* GM 



via Proximity Effect calculations by respective factors of roughly 1.9 and 2.3 ( Espey 1993 ). In 
addition, the lower limit to the density of baryonic matter in the Universe fib, is proportional 
to the square root of Jl GM ( [Rauch et al. 19971) . Therefore, over-estimating the strength of the 



background flux translates into an over-estimation of the lower limit of f2{, by factors of about 1.4 
and 1.5. 



We are grateful to Ray Weymann who generously reviewed this Letter and provided insightful 
comments and suggestions. We acknowledge helpful discussions with Robert Antonucci, Jill 
Bechtold, John Cocke, Mike Corbin, Rolf Kudritzki, Phil Pinto, Jennifer Scott, Matthias Steinmetz 
and Belinda Wilkes. And we thank the anonymous referee for fair and thoughtful suggestions 
that improved our Letter. This research was supported, in part, by National Science Foundation 
Grants AST 9529190 and 9803072. 
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Fig. 1. — Distributions of the broad emission line velocity shifts relative to the systemic frame 
(given by the [O III]A5007 emission line center). A negative velocity shift corresponds to a blueshift 
relative to systemic, while a positive shift represents a redshift. The velocity bins are ~ 1(a) for 
the broad H/3 shifts and ~ 2.5(a) for the UV rest-frame line shifts. The unweighted average of each 
blueshifted distribution is plotted as a solid dot, while the weighted mean with ±la uncertainty is 
shown for the redshifted broad H/3 distribution. The ionization potential Xion of each UV line is 
given, as well as its rest wavelength (in A). 
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Fig. 2. — The comparison between the distributions of broad H/3 velocity shifts for the nearby 
(z < 0.5, upper panel), lower luminosity sample of Boroson & Green (1992); and our ~ 100 
times more luminous, high redshift (z > 2.0, lower panel) sample. The broad component H/3 line 
centers, used to calculate the velocity shifts for each sample, were measured using the identical 
procedure (described in §2.2) and [O III] EW > 5A selection criteria. The mean rest-frame F-band 
luminosities of these two samples differ by roughly a factor of 100 (see M99). The data bin size 
is ~ 1(a) of the broad H/3 velocity shifts found in the more luminous sample. The mean of each 
distribution is plotted as a solid dot. 



